The  woth  re  pin-ice)  in  this  darmmnl  wrts  (itrlormcd  it  Lincoln  Laboratory,  a  <  enter 
tot  research  Of*r«n)  by  Massach  write  Institute  «i  Technology,  with  Jn  suppwt  of 
the  Department  of  the  Air  Force  under  Contract  FlMB-MK  -OWf 

Thia  report  may  b»  reproduced  to  rttlafv  nerd*  of  US.  Government  agencies 


The  view*  an  <  ai  m  contained  In  thlt  document  are  tbc'e  of  thr  contractor 
and  dw»M  *«  4  m*d  aa  aecttniily  tep  resenting  the  official  pvllcin,  riihn 

rapreaard  or  hr  nth-  •.  •  1  thr  United  State*  Cwtrwwwi. 


The  PvibUc  ‘M*Vv  h“»  »vi«rel  thit  tenon,  snrl  it  it 

releasable  to  the  Naiionai  Technical  Information  Service, 
where  ii  will  he  available  to  the  general  public,  including 
foreign  nationals. 


Tbi*  technical  report  has  been  reel  erred  and  ie  approved  lor  publication 
FOh  HF  f  OMMANDE.it 
A^Hw-r-se#7  ^  £****£& 

Rtyrrou d  l„  loiselle,  Lt.  Col.,  I-SAF 
Chief.  V..S4 >  t  Iim  c!ii  tdihorii>>ry  Ptojcct  Olluc 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 

OPTIMAL  SEARCH  THEORY  AND  OPTICAL 
ARTIFICIAL  SATELLITE  SEARCHES 


L.G.  TAFF 
Group  'M 


PROJECT  REPORT  ETS-64 

S  MAY  1982 


Approved  (or  public  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


- .. 'VsarWn  ,,, 

_■  *  *-Vi'  •'  _ _ 


ABSTRACT 


This  Report  discusses  the  along  orbit  search,  by  optical  means,  for  an 
artificial-  satell ite . . _ In  particular  the  attempt  is  made  to  couch  the  search 
In  the  existing  scenario  of  optimal  search  theory'.  This  cannot  ber  tione  forr - 
existing  and  envisaged  searches.  The  reasons  for  this  are  explored  and  some 
new  concepts  of  optimality  are  discussed  for  real  searches.  The  point  Is 
made  that  both  hardware  and  software  adjustments  would  be  necessary  In  order 
to  reconfigure  optical  searches  for  artificial  satellites  so  that  search 
theory  can  be  utilized. 
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I.  INTRODUCTION 

The  theory  of  search,  as  a  branch  of  operations  research,  had  its 
beginnings  during  World  War  II.  This  start  was  made  by  Bernard  Koopman  and 
his  colleagues  for  the  United  States  Navy  on  problems  related  to  anti¬ 
submarine  warfare.  Since  then  there  has  been  considerable  progress  in  the 
theoretical  development  of  searches  for  non-moving  targets.*  Some  progress 
Ras  Wen 'inicfe  "when ‘tHe  target  "is  allowed  to  move  in  a  conditionally  deter-'" 
ministlc  fashion.  Lawrence  D.  Stone  has  summarized  much  of  the  field  In  his 
book,  "Theory  of  Optimal  Search"  published  in  1975  by  Academic  Press.  This 
Report  is  an  introduction  to  this  field  of  mathematics.  I  closely  follow 
Stone's  notation.  This  is  to  enable  the  interested  reader  to  make  a  smooth 
transition  to  the  literature  should  he  care  to  pursue  the  subject  further. 

This  Report  Is  not  a  re-writing  of  Stone's  book.  Although  I.  follow 
him  in  basic  definitions  and  notation,  I've  chosen  to  illuminate  the  concepts 
within  the  framework  of  an  along  orbit  search  for  an  artificial  satellite  by 
optical  sensors.  As  the  Ground-Based  Electro-Optical  Deep  Space  Surveillance 
(GEODSS)  system  comes  on-line,  optical  searches  for  artificial  satellites 
will  become  more  frequent.  An  along  orbit  search  is  a  particularly  simple 
search  that  can  Illustrate  how  optical  searches  for  artificial  satellites  fit 

o 

Into  the  existing  theoretical  framework. 

Actually  such  searches  don't  fit  within  the  existing  framework.  A  minor 
problem,  for  rapid  searches  on  slow  moving  satellites,  is  the  fact  that  the 
target  Is  moving.  When  the  search  is  executed  slowly  or  the  satellite  Is 

*  I  deliberately  avoid  using  the  adjective  stationary  because  of  the  special 
meaning  it  has  in  the  artificial  satellite  context. 


moving  quickly,  this  becomes  a  major  problem.  The  principal  reason  that 
the  mathematical  superstructure  of  optimal  search  theory  is  superfluous  is 
that  the  optimization  problem  has  been  co-opted  by  the  basic  design  of 
existing  and  currently  contemplated  searches.  I  will  consider  some  of  the 
reasons  for  this.  This  Report  can  serve  as  a  stimulus  for  re-thinking  why 
we  search  the  way  we  do  and  whether  or  not  we  should  modify  It  to  obtain 
the  benefits  of  the  existing  theory. 

Given  the  existence  of  Stone's  book  I  won't  reference  any  other  litera¬ 
ture  herein.  All  results  specifically  pertinent  to  the  along-orbit  optical 
search  for  an  artificial  satellite  are  original. 
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II.  BASIC  IDEAS 


First  we  introduce  the  concept  of  a  discrete  search  space.  The  target 
is  In  one  of  a  set  J  of  cells  that  are  taken  to  be  a  (possibly  infinite)  sub¬ 
set  of  the  positive  Integers.  We  can  obtain  a  concrete  realization  of  this 
in  the  along  orbit  optical  search  problem  by  imagining  that  the  sensor  has  a 
field-of-view  of  9  radians  and  that  our  search  lies  along  the  artificial 
satellite's  orbit.  The  maximum  number  of  cells  is  2tt/9.  Geometrical  con¬ 
straints  would  limit  this  further.  Also  most  along  orbit  searches  are  based 
on  the  assumption  that  the  satellite  will  be  "early"  by  a  maximum  amount  and 
"late"  by  a  maximum  amount.  This  too  would  limit  the  maximum  number  of  cells. 

The  second  concept  is  that  of  an  a  priori  target  distribution.  By  this 
we  mean  the  probability  distribution  over  the  search  space  which  summarizes 
our  knowledge  of  where  the  target  was  likely  to  be  when  the  search  commenced. 
We  symbolize  the  target  distribution  by  p(j)  where 

P(j)e[0,l]Vjed  and  £  p(j)sl 
jed 

In  general  we  allow  the  target  distribution  to  be  defective.  In  the  along 
orbit  problem  we  might  have  to  allow  for  a  maneuver  of  the  satellite  which 
changed  its  orbital  plane.  If  our  search  space  consisted  of  only  a  few  cells, 
then  the  satellite  might  not  appear  in  any  of  them. 

When  we  search  a  cell  jeO  we  expend  a  certain  amount  of  effort  measured 
by  ze[0,°°).  The  cost  of  this  effort  is  given  by  the  cost  function  c(j,z). 

Cost  may  be  measured  in  time  or  money  while  effort  might  represent  area 
searched  or  the  duration  of  time  spent  in  a  cell.  In  the  along  orbit  problem 
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effort  and  cost  might  be  both  measured  in  time.  The  effort  would  be  simply 
the  amount  of  time  devoted  to  searching  cell  j.  The  cost  would  be  this  time 
plus  the  amount  of  time  necessary  to  move  the  telescope  to  cell  j  from  its 
last  location.  Such  a  cost  function  Is  realistic  and  does  not  fit  Into  the 
above  framework.  If  the  amount  of  time  spent  searching  in  a  cell  was  much 
longer  than  the  amount  of  time  required  to  move  to  a  cell,  then  we  might 
ignore  the  non-local  nature  of  the  cost  function.  We  could  for  instance,  use 
a  cost  function  of  the  form  c(j,z)  -  c  j  +  z  where  Cj  is  the  average  time  to 
move  to  cell  j  from  any  other  cell  In  J.  At  this  early  stage  we've  been 
forced  to  make  two  modifications  of  real  along  orbit  searches  in  an  attempt 
to  fit  Into  the  formalism;  real  satellites  are  moving  not  fixed  and  real 
cost  functions  are  non-local. 

When  we  expend  effort  z  in  cell  j  we  spend  cost  c(j,z).  In  return  for 
this  we  Increase  (or  at  least  not  decrease)  the  probability  of  detection. 

We  define  a  detection  function  b(j,z)  that  measures  this.  Specifically 
b(j,z)  is  the  conditional  probability  that  after  expending  effort  z  in  cell 
j  we  will  detect  the  target  given  that  It  is  in  cell  j.  Note  that  this 
definition  of  the  detection  function  assumes  that  the  conditional  probability 
of  detection  depends  only  on  the  total  amount  of  effort  expended  and  not  on 
the  way  that  the  effort  was  expended.  An  obvious  limitation  we  must  Impose 
on  the  cost  function  c(j,z)  is 

c(j,z)>0  v  jeO,  ze[0,») 


4 


Clearly  we  must  demand  that  the  detection  function  satisfy 
b( j ,z)  e[0,l]  V jeJ,  zc[0,°°) 

There  is  a  special  class  of  detection  functions  known  as  regular  detection 
functions  that  are  important  in  the  theoretical  framework.  The  detection 
function  b(j,z)  is  regular  if 

(i)  b( j ,0)  =  0  and 

(11)  db(j,z)  is  continuous,  positive,  and  strictly  decreasing  \/je.J 
<?z 

The  first  requirement  just  says  that  we  must  expend  effort  If  we  are  to  find 

the  target.  The  second  requirement  implies  that  b(j,z)  is  strictly  concave 

and  that  the  rate  of  return  function  p(j,z)  =  p(j)db( j ,z)  diminishes  with 

6Z 

increasing  effort.  The  rate  of  return  function  tells  us  which  cell  yields 
the  largest  increase  in  probability  of  detection  for  a  given,  small,  incre¬ 
ment  of  effort.  A  search  that  always  places  the  next  increment  of  effort  in 
the  cell  with  the  highest  rate  of  return  Is  called  locally  optimal . 

The  detection  function  b(j,z)  Is  something  that  we  could  model  for  a 
particular  optical  sensor.  Presumaoiy  it  would  require  only  minor  approxi¬ 
mations  (or  perhaps  no  approximations)  for  such  a  detection  function  to  be 
regular.  As  far  as  I  am  aware  this  has  not  been  done.  It  would  turn  out  to 
be  superfluous  for  presently  conceived  searches  though  because  current  hard¬ 
ware  constraints  are  such  that  one  would  not  change  the  amount  of  effort 
allocated  to  a  cell  during  a  search.  This  ab  initio  fixing  of  the  effort 
allocated  to  a  cell  (in  fact  to  be  the  same  constant  amount  for  each  cell) 


co-opts  the  optimization  problem  completely.  We'll  see  this  more  clearly 

soon.  Moreover,  even  if  this  were  not  the  case,  many  types  of  artificial 
satellite  searches  have  designs  which  prohibit  the  design  of  an  optimal 
search  in  the  sense  discussed  below.  Searches  that  define  J  or  searches  that 
are  designed  to  be  leakproof  can  remove  the  essential  degree  of  freedom 
necessary  for  the  successful,  meaningful  application  of  search  theory. 

So  far  we  have  introduced  the  search  space  0  of  cells,  the  a  priori 
probability  distribution  over  the  search  space  p(j),  the  cost  of  expending 
effort  z  in  cell  j,  c(j,z),  and  the  conditional  probability  of  detecting  the 
target  in  cell  j  after  expending  effort  z  there,  b(j,z).  Finally  we  intro¬ 
duce  the  concept  of  an  al location  of  effort  over  J,  f(j).  If  we  define  the 
set  F(J)  by 

F  =  F( J)  =  Set  of  all  non-negative  functions  f  defined  on  J 

then  we  can  formulate  the  basic  search  problem. 

Suppose  we  have  an  allocation  of  effort  f(j)  for  each  cell  in  the  search 
space.  Then  the  total  cost  C[f]  of  this  allocation  Is 

C[f]  =  Z  c(j,f(j)) 
jeJ 

The  total  probability  of  detecting  the  target  with  this  allocation  is  P[f], 


P[f]  =  Z  P( j )b{ j  ,f ( j ) ) 
jeJ 


Suppose  that  the  maximum  cost  is  K.  The  basic  search  problem  is  to  find 
an  allocation  of  effort  f*eF  such  that 

C[f*]<  K 
and 

P[f*]  =  max  | P[f 3 :  feF  and  C[f]  s  <J 

If  such  an  allocation  exists  then  it  Is  said  to  be  optimal  for  cost  K. 

Now  we  can  clearly  see  the  critical  element  of  the  above  mentioned 
difficulty.  We  ignore  the  fact  that  the  theory  is  concerned  with  fixed 
targets.  We  ignore  the  fact  that  the  theory  is  capable  of  dealing  only 
with  local  cost  functions.  We  cannot,  however,  ignore  the  fact  that  opera¬ 
tionally  we  ab  initio  determine  that  f(j)  is  a  constant  such  that  C[f]  =  K. 
This  stricture  removes  the  free  element  of  the  optimization  problem.  To  fit 
optical  artificial  satellite  searches  within  the  theoretical  framework  we  can 
not  ab  Initio  declare  that  f(j)  is  a  constant.  In  order  to  free  the  alloca¬ 
tion  of  effort  requires  a  re-design  of  some  hardware  and  software  as  well 
as  careful  examination  of  the  current  demands  on  artificial  satellite  searches 
(eg  leakproofness). 

After  we  introduce  one  more  concept  we  can  turn  to  the  results  of 
optimal  search  theory.  The  final  concept  is  that  of  a  search  plan.  A  search 
plan  4>(j,t)  Is  a  function  defined  on  jeJ,  tc[0,“)  which  tells  us  how  much 
effort  has  been  expended  in  cell  j  by  time  t.  If  M(t)  gives  the  total  effort 
available  by  time  t  then 

£  <t>( j  *t)  =  M(t) 

jeJ 
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Let  $>(M)  be  the  . '  ass  of  search  plans  satisfying  the  above.  Then  we  say  that 
a  search  plan  is  uniformly  optimal  within  4>(M)  if 

P[4>*(j.t)]  =  max  |p[4>(J*t)3: <t>e4> ( M ) J  v  t>0 

In  other  words  such  a  search  plan  maximizes  the  probability  of  detection  at 
every  instant.  Under  certain  conditions  one  can  show  that  the  locally  optimal 
plan,  the  uniformly  optimal  plan,  and  the  plan  that  minimizes  the  mean  time 
to  detection  are  identical. 


III.  THEORY 

Above  the  basic  concepts  of  search  theory  in  a  discrete  search  space 
have  been  Introduced.  In  this  section  we  shall  see  how  optimal  search  plans 
can  be  found  and  computed.  No  proofs  are  provided;  they  are  in  Stone's  book. 
Also,  as  I  feel  that  it's  more  important  to  convey  the  sense  of  what  can  be 
demonstrated,  rather  than  the  details,  I've  taken  some  liberties  in  relaxing 
the  wording  of  some  of  the  results.  Again  see  Stone's  book  for  details. 
definition:  An  allocation  of  effort  f*eF(J)  is  optimal  for  cost  K  if 

C[f*]<K  and  P[f*]  =  max  |p[f]:feF(J)  and  C[f 3 S k| 

Theorem:  Suppose  there  exists  a  Xe[Q,»)  and  an  allocation  f*eF(J)  such 
that  C[f*]<-  and 

p[f£]  -  XC[f*J>P[f]  -  XC[f]  for  feF(J)  3C[f]<*> 

then 


P[f*]  =  max  |p[f]:feF(J)  and  C[f]SC[f*]| 

This  theorem  tells  us  when  an  allocation  fj  is  optimal  for  a  given  cost 
C[f*].  The  Introduction  of  the  Lagrange  multiplier  X  converts  a  constrained 
optimization  problem  for  the  functional  P[f]  into  an  unconstrained  optimization 
problem  for  P[f]-XC[f].  This  device  suggests  the  introduction  of  the  point- 
wise  Lagrangian 

£U,X,z)  *  p( j )b( j ,z )-Xc(j  ,z)  for  V  jeJ;  X,ze[0,°°) 
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!  1 

i  ; 

s 

Note  that 

E-C(j,X,f)  -P[f]-XC[f]  I 

jeJ  | 

Theorem:  Suppose  there  exists  a  Xe[0,«)  and  an  allocation  f£eF(J)  such  that  j  | 

C[f£]<*  and  !  j 

>  \ 

£(j,X,f*(J))  «  max  |£{j ,X.z):ze[0,«)|  VjcJ  i  ] 

|  v; 

then 

P[f*]  »  max  |p&f]:feF(J)  and  C[f]<C[f£]|  j  ! 

I 

} 

j . 

Note  that  in  this  simple  case  once  j  and  X  are  fixed  t(j,X,z)  Is  a  function  j 

of  one  variable  and  the  apparatus  of  ordinary  differential  calculus  may  be  j  , 

|  ; 

used  to  find  Its  maximum.  The  above  theorems  give  sufficient  conditions  i 

j 

that  a  particular  allocation  satisfies  the  constrained  optimization  problem.  ! 

I 

To  go  further  we  must  assume  a  bit  more  about  the  detection  function  and  the  j; 

cost  function.  Moreover  we  say  that  whenever  (X,f*)  are  such  that  f*eF(J) 
and  Xe[0,«)  and 

£(j »X,f*(j ) )  *  max  |£(j,X,z):ze[0,°°]  j  VjeJ 

I 

i 

that  (X,f*)  maximizes  the  pointwise  Lagrangian. 

Theorem:  Let  b(j,z)  be  a  concave  function  of  z  and  c(j,z)  be  a  convex 

function  of  z.  Let  f*eF(J)  and  C[f*]e(Q,K).  Then  a  necessary  and  sufficient 

condition  for  f*  to  be  optimal  for  cost  C[f*]  is  that  there  exists  a  Xe[0,°°) 

such  that  ( X , f * )  maximizes  the  pointwise  Lagrangian.  I 
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Corollary:  Let  db(j,z)/dz  be  a  decreasing  and  continuous  function  of  z  and 
let  dc(j,z)/dz  be  an  increasing  and  continuous  function  of  z  vjeO.  If 
f*eF(J)  Is  optimal  for  cost  C[f*]e(0,K)  then  there  exists  a  X>0  such  that 

*  0  If  f*(j)  >  0 

P(j)  3b(J,f*(J))  -  X  3G(J.f*(J)) 

dz  dz  <  0  if  f*(j)  *  0 

To  complete  the  elementary  part  of  the  theory  we  need  the  Inverse  of 
the  rate  of  return  function  p(j,z)  *  p(j)  <3b(j  ,z)/<3z,  viz. 

p-i(j  -  Inverse  of  p(j,z)  evaluated  at  zax  for  xe(0 ,p(j  ,0)] 
0  for  x  >  p(j ,0) 


We  also  need  the  function 


U(x)  »  £  p_1(J,x) 
jeJ 


and  Its  inverse  U  • 

Theorem:  If  c(j,z)=z  \j ze[0,°°)  and  jeJ  and  b(j,z)  is  a  regular  detection 
function  then  for  a  fixed  cost  K>Q  the  allocation 


fj  -  P'1  (j,X)  VjeJ 


where  X  *  U~ 1  ( K)  is  optimal  for  cost  K  and  C[f£]  ■  K. 
Theorem:  Under  the  above  conditions  the  search  plan 


4*  (J.t)  -  P'1  (j  «U"1  (M(t))j  VjcO,  t>0 


is  uniformly  optimal  for  cumulative  effort  M(t)  in 
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Finally  define  the  mean  time  to  find  the  target  when  using  search  plan  y  (<*>). 
We  can  write 

p(4>)  3  Jq  (l-PU(j.t)])dt 

Theorem:  Let  the  search  plan  41*  be  uniformly  optimal  In  4>(M).  Then 

C  4>)  • 

The  above  theorems  completely  solve  the  problem  of  constructing  optimal 
search  plans  for  regular  detection  functions  and  cost  functions  proportional 
to  the  effort  expended.  In  order  to  give  some  life  to  theory,  consider  the 
along  orbit  search  again.  Let  there  be  J  cells,  let  the  target  distribution 
be  p(j),  and  suppose  that  the  detection  function  b(j,z)  =  l-exp(-ajZ) ,  a^O. 
Note  that  this  detection  function  Is  regular.  Fix  the  total  cost  at  K  and 
let  the  cost  function  be  c(j,z)=z.  We  will  construct  an  algorithm  for 
finding  the  optimal  allocation  of  effort  f{d)32j. 

The  total  probability  of  detection  for  allocation  f  is 

J  J 

p[f]  =  I  P(j)b(j,f(j))  -  E  p(j)0-exp(a  .z  .)) 
j-1  J-l  J  J 

The  total  cost  for  this  allocation  Is 

J  J 

C[f]  ■  £  C(j,f(j))  .£  z j 
j»l  j-l  J 

We  seek  to  maximize  P[f]  subject  to  the  constraint  C[f]*K.  We  do  this  by 
introducing  the  Lagrange  multiplier  X  and  minimizing  P[f]-AC[f]  with  respect 
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to  {z  . } .  We  find  that  <?(P[f]-XC[f  ])/<3z  .=0  implies 

J 

p(j)ajexp(-ajZj)  =  X  Vje[l,J] 

Hence 

Zj  *  O/aj)  ir\  (p(j)ajA) 

Notice  that  the  implicit  equation  for  Zj  is  just  p(j,Zj)  =  X  and  that  the 
solution  is  just  Zj  =  p“ 1 ( j ,X ) .  Of  course,  we  can  only  consider  non-negative 
effort  so  if  p( j)aj/X<l,  Zj=0.  We  determine  X  from  the  total  cost  con¬ 
straint, 


J  J 

Lzj  “  K=  £  0/a.j)  (P^)aj/X)  a 

j-1  J-l 

where  only  those  terms  are  included  whose  argument  of  the  natural  logarithm 
are  21.  This  completes  the  solution  of  the  problem.  To  see  even  more 
clearly  how  the  search  evolves  suppose  that  the  cells  are  labeled  such  that 

a1  p(  1 )  2a2p(2) . >  a0p(J) 

If  the  total  cost  K  Is  small  only  cell  #1  will  be  searched  for  p(j)a^  will 
be  less  than  X  for  j>l.  As  the  total  cost  increases  additional  cells  will 
be  searched  with  the  effort  partitioned  between  them  according  to  the  above 
rules. 
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IV.  ADDITIONAL  TOPICS 

The  above  theoretical  development  concerned  optimal  search  plans  over 
a  discrete  search  space.  The  amount  of  effort  expended  was  infinitely 
divisible  though.  A  more  realistic  approach  would  be  to  quantify  the  amount 
of  effort  expended  in  each  cell.  All  expenditures  of  effort  would  then  occur 
in  discrete  multiples  of  the  minimum.  In  optical  searches  the  minimum  amount 
of  effort  is  the  time  to  form  an  Image  with  the  sensor.  Clearly  it  makes  no 
sense  to  allocate  an  amount  of  effort  (measured  in  time)  that  is  less  than 
the  single  image  integration  time. 

Such  searches  are  known  as  search  with  discrete  effort.  A  logical 
measure  of  the  cost  is  the  amount  of  time  needed  to  perform  a  single  look 
(form  an  Image)  in  a  cell.  The  quantity  subject  to  variation  in  the  opt1mi: 
zation  problem  is  the  distribution  of  the  number  of  looks  per  cell.  This  is 
analogous  to  the  allocation  of  effort  discussed  above.  One  can  also  carry 
over  the  notion  of  a  locally  optimal  search.  Such  a  search  plan  looks  In 
that  cell  that  yields  the  highest  value  for  the  quantity  (increment  in  pro¬ 
bability)/  (increment  in  cost).  If  this  ratio  decreases  with  the  number  of 
looks  (in  every  cell  of  the  search  space),  then  the  locally  optimal  plan 
minimizes  the  mean  time  to  find  the  target.  If  In  addition  we  do  measure 
cost  by  the  number  of  looks  per  cell,  then  the  locally  optimal  plan  is  also 
uniformly  optimal,  i.e.  It  maximizes  the  probability  of  detection  for  any 
number  of  looks  >0. 

Other  topics  treated  by  the  theory  include  whereabouts  searches, 
optimal  search  and  stop,  search  in  the  presence  of  false  targets,  the  approx¬ 
imation  of  optimal  search  plans,  some  small  steps  In  solving  the  conditionally 
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deterministic  target  motion  problem,  and  Markovian  target  motion. 
Interested  reader  Is  referred  to  Stone's  book. 


The 
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V.  THE  ALONG  ORBIT  SEARCH 


Within  the  framework  of  currently  conceived  searches  there  would  appear 
to  be  little  left  to  discuss.  One  can  still  pose  interesting  problems  for 
along  orbit  searches  depending  on  one's  concept  of  optimality.  In  this 
Section  I'll  look  at  some  of  the  simple  problems  an  along  orbit  search 
presents. 

A.  Search  Scenario 

The  search  Is  cv;ns trained  in  the  following  fashion.  The  search 
space  is  a  set  of  2N+1 ,  N 2t0  (an  Integer)  fleld-of-view  0.  The  fields  are 
labeled  n«-N,  -N+l,  ....  N.  All  searches  commence  at  n*0,  the  nominal 
position  for  the  artificial  satellite.  The  2N+1  fields  lie  along  the 
satellite's  orbital  plane.  The  target  distribution  p(n)  is  specified.  The 
restrictions 

1  £p(0)  t  p(n)  SO  V  ne[-N,N] 

might  be  Imposed  with  little  loss  of  generality.  In  addition  one  might 
Impose  a  symmetry  constraint  p(n)  »  p(-n) y  ne[-N,N]. 

Each  cell  Is  searched  with  the  same  amount  of  effort.  Moreover 
we  assume  that  the  detection  function  Is  homogeneous,  b(n,z)  *  b(m,z) 

\/  n,me[-N,N].  We  define  a  search  plan  as  a  set  of  2N+1  Integers  n0,n1,n2  ... 
n2fj  drawn  from  [-N.N]  subject  to  the  constraints  that  nQ  *  0  and  that  there 
is  no  repetition.  There  are  (2N)1  different  search  plans  of  which  half  are 
the  reflection  of  the  other  half.  The  cost  function  we  develop  below  after 
we  construct  a  model  for  the  telescope's  motion.  Then  we  can  pose  several 
questions:  1)  Which  search  plan  takes  the  least  time  to  complete  an 
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examination  of  all  2N+1  cells?,  2)  Which  search  plan  has  the  highest  average 
probability  of  finding  the  target  as  a  function  of  time?  3)  Which  search 
plan  has  the  highest  aggregate  probability  of  detection  for  all  times  during 
the  search?  (The  aggregate  probability  of  detection  is 

P(t)  *£p(n) 
neM(t) 

where  M(t)  Is  the  set  of  cells  searched  by  time  t.),  and  4)  Do  any  of  these 
matter  In  real  world  along  orbit  optical  artificial  satellite  searches? 

B.  A  Model  for  the  Telescope 

In  order  to  define  a  realistic  cost  function  we  need  a  model  for 
the  telescope  motion.  For  this  purpose  I  assume  that  when  the  telescope 
starts  from  rest  that  it  Is  capable  of  a  maximum  constant  acceleration  a 
for  a  maximum  time  t.  Hence  the  maximum  angular  speed  of  the  telescope  is 
n  =  at.  The  telescope  can  move  at  the  rate  of  D  for  an  arbitrary  length  of 
time.  When  the  telescope  decelerates,  it  does  so  at  a  constant  deceleration 
6  ( <3 > 0 )  until  It  comes  to  rest  or  a  specified  angular  speed  w, |w|£fi. 

Suppose  that  the  problem  Is  to  move  the  telescope,  initially  at 
rest  at  cj)  ^  to  some  other  position  where  It  will  again  be  at  rest. 

This  process  can  occur  in  a  maximum  of  three  phases.  During  phase  I  the 
telescope  accelerates  at  the  rate  of  a  for  a  total  time  t  <x.  During  phase  II 
the  telescope  moves  at  the  constant  angular  speed  w  *at.  for  a  time  t. 

During  phase  III  the  telescope  decelerates  at  the  rate  6  until  it  reaches  a 
stop  at  4>  This  takes  an  addition  time  t^.  The  total  time  of  the  move 
Is  T  *  ta+t-+td.  Question:  Wnat  combination  of  ta,  tc  and  t^  minimizes  T? 


17 


To  solve  this  problem,  we  need  equations  of  motion  for  the  telescope. 
They  are 

Phase  I:  te[Q,ta] 

<j>  «a  t2/z  +  <J>.| »  uj  =  a  t 

Phase  II:  t-ta  e[0,tc]  4>a  !  «  t|i  +  ^  ,  uf  ata 
<!>  *  ^  (t-ta)  +  4>a  .  w  »  wa 

Phase  III:  t-(ta+tc)  e[0,td]  $c  -  <^a  tc+  4>a  .  “c  =  «a 
<0  *  -6[t-(ta+tc)]2A  +  uc[t-(ta+tc)]  +  *c  . 
io  ■  -5  [t-(ta+tc)3  +  uc 

Set  t  -  t  +t  +trf  In  the  Phase  III  equations  and  Insist  that  $  «  <t>f  .  w  *  0. 

d  C  Q 

One  finds 

td  -  a  tfl/S  ,  tc  (l+a/6) 

where  A<fr ■  <J> f  -  ^ .  The  total  traverse  time  T  Is 
T  -  ta  (l+a/6)  +#■ 

Considered  as  a  function  of  tg,  T  has  a  single  minimum  when 

*2  '  2A<t> 

*a  “  a(l+a/,$  ) 

which  Implies  that  the  coast  time,  tc  Is  zero.  The  total  time  to  move  an 


angular  distance  A<j>  is 


The  above  assumes  that  t  <  t.  Should  the  above  value  of  t  be 

a  a 

greater  than  or  equal  to  x,  then  the  single  degree  of  freedom  is  removed  and 


Only  for  an  extremely  fast  telescope  or  a  very  long  along  orbit  search 


would  t  exceed  x. 

a 


C.  Optimal  Searches 

The  cost  function  we'll  use  is  the  time  required  to  move  from  the 
last  cell  examined  to  the  cell  of  Interest  plus  the  time  spent  examining  the 
cell.  The  time  per  cell  (=»  the  effort  expended)  is  a  constant  equal  to 
tiook*  Since  the  center  of  cell  ne[-N,N]  is  at  41=  ne,  for  tfl<x  the  cost 
function  for  the  K'th  cell  is 


'  nK-nK-l  l]  /2  +  ^ok 


Consider  first  a  complete  search  (1,e.  each  cell  is  examined).  Then  the 
total  time  to  complete  the  search  is 


(2N+1) tiook 


j2(cth$)8  j 


lnl"n0l  +  |iYV  +  .  . 


+  I^N^N-l 


Clearly  when  determining  an  optimal  search  only  the  sum  of  the  square  roots 
are  Important. 


Conjecture:  The  complete  search  plan  that  has  the  shortest  time  to  complete 
is  the  following  (or  its  reflection):  Pick  a  direction  and  move  in  it,  one 
cell  at  a  time,  until  the  last  cell  on  that  side  is  reached.  Next  jump,  in 
the  other  direction,  one  cell  past  the  start  and  repeat  the  (new)  unidlrec- 
tional  one  cell  at  a  time  traverse.  The  time  to  complete  is  2N+(N+1 )  /2-l . 

The  along  orbit  searches  actually  used  is  an  alternating  one,  eg.  nQ=0, 

n^3+l,  n23-l,  n3*+2,  n^=  -2,  etc.  The  time  to  complete  Is 

2N 

L  n*  ^N3/2 

n*l 

and  therefore  very  long.  This  search  does  build  up  aggregate  probability 
quickly  though,  especially  for  a  sharply  peaked,  unimodal ,  symmetric  distri¬ 
bution.  If  the  arc  of  the  orbit  is  long,  then  this  probably  isn't  optimal 
in  the  sense  of  maximum  aggregate  probability.  The  reason  is  that  since 
p(n)  falls  off  rapidly  with  n  but  the  time  to  move  Increases  as  the  square 
root  of  n,  for  large  enough  of  n  it  will  be  better  to  do  cell  n+1  after 

cell  n  and  then  jump  to  cell  -n,  do  cell  -(n+1),  then  jump  to  cell  n+2,  etc. 

Clearly  one  needs  real  numbers  for  a,6,6,  and  t^Qok  to  decide  the  question. 
The  same  is  true  for  the  highest  average  probability  searches. 

D.  Does  it  matter? 

Stone  references,  but  does  not  deal  with,  the  subset  of  the  search 
literature  concerned  with  searches  along  a  line.  I've  briefly  looked  at  it 
and  It  appears  to  be  Irrelevant.  Finally,  unless  one  has  an  extremely  slow 
telescope  or  contemplates  very  long  (in  the  sense  of  arc)  along  orbit 
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searches,  it  is  doubtful  that  considerations  of  optimality  really  matter. 
Note  though  that  implicit  in  current  along  orbit  searches  is  the  assumption 
that  p(0)  n,  Vz.  P (i  1 )  ^  Yu,  or  the' satellite  maneuvered. 
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